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Kinematic distributions in the decays of the newly discovered resonance to four leptons are a 
powerful probe of the tensor structure of its couplings to electroweak gauge bosons. We present 
analytic calculations for both signal and background of the fully differential cross section for the 
'Golden Channel' e + e~ fi + fi~ final state. We include all interference effects between intermediate 
gauge bosons and allow them to be on- or off-shell. For the signal we compute the fully differential 
decay width for general scalar couplings to ZZ, 77, and Z7. For the background we compute the 
leading order fully differential cross section for qq annihilation into Z and 7 gauge bosons, including 
the contribution from the resonant Z — > 2e2/i process. We also present singly and doubly differential 
projections and study the interference effects on the differential spectra. These expressions can be 
used in a variety of ways to uncover the nature of the newly discovered resonance or any new scalars 
decaying to neutral gauge bosons which might be discovered in the future. 



I. INTRODUCTION 

With the recent discovery of a new resonance at the 
LHC [TJ [2] the focus now shifts to the determination of 
its detailed properties including its spin, CP, and elec- 
troweak (EW) quantum numbers. It has been shown in 
recent studies 3 8 and also emphasized for quite some 
time [9 - 13J , that the decay to four charged leptons is 
a powerful channel in accomplishing this goal. Because 
of the experimental precision with which this channel is 
measured, it offers one of the few opportunities to use an- 
alytic methods to analyze the data. We thus seek to ex- 
tend previous analytic calculations of both the signal and 
the standard model (SM) background and present com- 
pletely general, leading order (LO) fully differential cross 
sections for the 2e2/i final state mediated by intermediate 
Z and 7 gauge bosons. In addition to performing discov- 
ery/exclusion analysis and signal hypothesis testing one 
could, with enough data, experimentally determine all 
possible couplings of a spin-0 scalar to pairs of neutral 
electroweak (EW) gauge bosons in one multi-parameter 
fit using these expressions. 

Analytic expressions are ideal for use in the matrix 
element method (MEM) taking full advantage of all of 
the kinematic information in the event. One can then 
use the fully differential cross section to build a likeli- 
hood function [14j [15] to be used as a discriminant. For 
a recent study of the golden channel comparing existing 
leading order MEM-based approaches and software [Po] . 
along with providing code which calculates kinematic dis- 
criminants based on the Madgraph [17] matrix element 
squared see [IS]. We view this 'analytic approach' as 
equivalent and complementary to these other approaches. 
These analytic expressions also allow for more flexibility 
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in performing multidimensional fits to determine cou- 
pling values which will be useful when performing pa- 
rameter extraction. Our parametrization allows for easy 
implementation of various hypothesis tests as well as the 
addition of NLO effects which can also be implemented 
into an MEM framework [T91I20] . 

For the signal we compute the fully differential decay 
width for the process ip — > ZZ + Z7 + 77 — > 2e2[i where 
ip is a spin-0 scalar. We allow for the most general CP 
odd/even mixtures and include all interference between 
intermediate vector bosons. While these expressions are 
applicable to the newly discovered boson at 125 GeV, 
they are also applicable for any new scalar decaying to 
neutral gauge bosons. This allows one to consider a va- 
riety of hypotheses which can be tested against one an- 
other. It should be emphasized however that for opti- 
mal performance, even when testing between two dif- 
ferent signal CP and spin hypothesis, one should also 
include the background in the discriminant since in any 
given sample it is not known with full certainty which are 
background and which are signal events. Thus we seek to 
provide both signal and background distributions which 
can be used together to build a complete likelihood. 

For the background we compute the fully differen- 
tial cross section for the qq 2e2[i process. This in- 
cludes the contributions from all the intermediate vector 
bosons through both t-channcl pair production and the 
singly resonant four-lepton production s-channel process 
qq — > Z — > 2e2/.i. We include all interference effects be- 
tween the intermediate vector bosons as well the interfer- 
ence between the s-channel and t-channel diagrams which 
can affect the differential distributions as we will see be- 
low. Also, unlike the analytic calculations in [5T] of 
the golden channel background differential cross section, 
these expressions are valid for a much larger energy range 
for the four lepton invariant mass as well as the invariant 
mass of each lepton pair. In particular, since these also 
include the 77 contribution one can probe lower values 
in the differential mass spectrums, which as we will see is 
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a highly discriminating region. The intermediate vector 
bosons are allowed to be on or off-shell and in what fol- 
lows we do not distinguish between the two. We do not 
discuss the 4e and 4/i final states explicitly, but in some 
kinematic regimes the interference effects between iden- 
tical final state particles can be sizable [18] . We leave an 
inclusion of these final states to future work. 

Although other channels can also probe the tensor cou- 
plings of a resonance to neutral gauge bosons, the golden 
channel, with a four body final state has the advantage 
of extra kinematic variables, such as the azimuthal angle 
between lepton decay planes. This variable would be un- 
available for example in the 77 final state. In addition to 
offering more kinematic observables, the golden channel 
offers the unique opportunity to test all of the possible 
tensor couplings including any potential interference ef- 
fects between the different operators in one direct (and 
very precise) fully correlated measurement without any 
recourse to theoretical input (other than the production 
cross section of course) . This allows for stringent tests of 
the SM to be performed and perhaps allow us to uncover 
new physics which may be hiding in subtle effects within 
the golden channel. 

In addition to presenting the calculation of the fully 
differential cross sections we examine various singly and 
doubly differential distributions and elucidate the subtle 
interference effects between the different contributions to 
the signal and background. Of course a proper treatment 
of the golden channel requires careful study of detector 
resolution and acceptance effects, but we leave that to 
ongoing analyses. 

The organization of this paper is as follows: in Sec. [TT] 
we briefly review the kinematics of the four lepton final 



state. In Sec.ylljwe describe the calculation of the signal 
fully differential cross section and examine the differential 
mass spectra for a variety of signal hypotheses. In Sec. |IV| 
we describe the calculation of the background fully dif- 
ferential cross section and examine how the kinematic 
variables are affected by NLO and pdf effects before con- 
cluding in Sec. |Vj We also present in the Appendix a 
pair of expressions for the signal and background dou- 
bly differential mass spectrums and also show plots for a 
multitude of singly and doubly differential spectra. 



II. FOUR LEPTON EVENTS 



to charged leptons one of which radiates a Z/7, which 
again decays to charged leptons (see Fig. [4} . The first 
lepton pair which radiates the second vector boson does 
not reconstruct to the Z boson four momentum, which 
in this case is also equal to the invariant mass of the At 
system. The kinematics remain unchanged, but now we 
must interpret the angles defined in the lepton pair rest 
frame with respect to the direction of momentum of the 
lepton pair system as opposed to that of one of the gauge 
bosons. Thus, we have the following more general inter- 
pretations for the kinematic variables defined in the At 
rest frame; 

• Mx 2 - The invariant mass of the two lepton pair 
systems. 

• 9 - The 'production angle' between the momentum 
vectors of the lepton pair which reconstructs to M± 
and the total At system momentum. 

• $1.2 ~~ Polar angle of the momentum vectors of 
e~ , /i - in the lepton pair rest frame. 

• $1 - The angle between the plane formed by the 
M\ lepton pair and the 'production plane' formed 
out of the momenta of the incoming partons and 
the momenta of the two lepton pair systems. 

• $ - The angle between the decay planes of the final 
state lepton pairs in the rest frame of the At system. 

These variables are all independent subject to the con- 
straint (Mi + M2) < y/s where s is the invariant mass 
squared of the 41 system. We have also ignored the irrel- 
evant azimuthal production angle. 

In the case of the signal events one can replace 'lepton 
pair' momentum with Z or 7 momentum since in those 
cases, both lepton pairs do indeed decay from a vector 
boson and the intuition follows that found in Fig [l] The 
same can be said for background events which proceed 
through t-channel pair production. In these cases, the an- 
gle $1 defines the azimuthal angle between the di-boson 
production plane and the plane formed by the lepton pair 
which reconstructs to M 1 and G is the vector boson pro- 
duction angle. Other than this more subtle interpretation 
of the various kinematic variable however, in practice the 
definitions of these variables are left unchanged from the 
definitions found in [3J which we follow from here on. 



The kinematics of four lepton (A£) events are described 
in detail in many places in the literature and here we 
use the convention found in [3J. We comment on the 
kinematics briefly and point out that in the case of the 
background the physical interpretation of the kinematic 
variables is not as intuitive as in the case of previous stud- 
ies which only considered the t-channel ZZ contribution. 
Now since we include the contribution from resonant four 
lepton production, the lepton pairs do not necessarily re- 
construct to a physical particle. In this case, resonant 
production of a Z (or possibly 7) is followed by decay 



III. SIGNAL 

In this section we present the calculation of the signal 
fully differential cross section and examine the differen- 
tial mass spectra for several signal hypothesis. We take 
our signal to be a general spin-0 scalar and consider all 
possible couplings to any combination of Z and 7 pairs 
allowing for mixtures of both CP even and odd interac- 
tions. Previous studies have analytically computed the 
ZZ [3J 3] contribution to the golden channel, but as far 
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FIG. 1. Definition of angles in the four lepton CM frame X. 

as we are aware, none consider the contributions from the 
Zj and 77 intermediate states. There are also interfer- 
ence effects between the intermediate state which are not 
present when 7 is not allowed to decay. As we will see, 
these effects can manifest themselves in the kinematic 
distributions. Of course for a SM Higgs, the Zj and 77 
contributions to the golden channel arc expected to be 
small, but this need not be true for a general scalar or 
if the discovered resonance turns out to have enhanced 
couplings to Z~f or to 77. How large these effects are 
once one takes into account detector and acceptance ef- 
fects deserves careful study, but we leave this for ongoing 
work. 

The most general couplings of a spinless particle to 
two gauge bosons with four momenta k\ and k 2 can be 
expressed as, 

iT% = v- x (A Uj m%gV + A 2ij (k 1 ■ k 2 g^ - k%k$) 

+A 3ij e IJ . l , o ,0k™k 2 3 ) (I) 

where ij = ZZ, Z7, or 77. The ^1,2,3 are dimensionless 
arbitrary complex form factors and v is the Higgs vacuum 
expectation value (vev), which we have chosen as our 
overall normalization. For the case of a scalar coupling 
to Zj or 77 electromagnetic gauge invariance requires 
A\ = 0, while for ZZ it can be generated at tree level 
as in the SM or by higher dimensional operators. We 
have chosen to write the vertex in this form to make 
the connection with operators in the Lagrangian which 
may generate them more transparent. For example the 
following list of operators may generate a coupling as in 
Eq.Q, 

C ~ ^(g h m%Z>*Z M + gzZ^Z^ + g z Z^Z^ 
+ g Zl F^Z^+g Zl F^Z^ 
+ g 1 F^F liV + ~g 1 F^F llv + ..)j (2) 

where Z^ is the Z field while = d^v—dyV^ the usual 



bosonic field strengths. The dual field strengths are de- 
fined as V^v = \(. llvpa V pa and the ... is for operators of 
dimension higher than five. For a given model many of 
these are of course zero. If <p is the Standard Model Higgs, 
then gh = i, while gz, gz-y and g 77 are 7^ 0, but loop in- 
duced and small. Detailed studies of the ZZ contribution 
to the golden channel mediated through the operators 
with coefficients gh, gz were conducted in [31 El 123 ■ T ne 
operators corresponding to gz-y were studied in [5] for the 
golden channel final state and in [23] for the £ + £~j final 
state and both were shown to be useful discriminators. 

Other recent studies of these operators, though not 
only through the golden channel final state, have also 
been done. The pseudo scalar couplings gz, gz-y, g-y were 
studied recently in the context of the newly discovered 
resonance in |24j where it was shown that a purely CP 
odd scalar is disfavored as the new resonance. The anal- 
ysis of [25] shows that with a fit of the 77, ZZ* , and 
WW* rates, as well as the absence of a large anomaly in 
continuum Z7, that the scenario of the four lepton de- 
cays being due to gz or gz-y is strongly disfavored. While 
these statements contain few assumptions, they are still 
model dependent and should be confirmed by direct mea- 
surements. 

Even if the newly discovered resonance appears to be 
'SM like', it is still possible that it can have contributions 
to the 2e2/i channel coming from operators other than gh 
which are slightly enhanced relative to the SM prediction. 
Here we are motivated by asking what information can 
be extracted from this channel with out any a-priori ref- 
erence to other measurements or theoretical input. In 
addition, there still exists the possibility that another 
scalar resonance will be discovered which can also decay 
to EW gauge boson pairs. In this case it may have com- 
parable contributions from the various operators. Thus 
we allow for all operators in Eq. ^ to contribute simul- 
taneously including all interference effects between the 
ZZ, Zj, and 77 intermediate states. Because the vertex 
in terms of arbitrary complex form factors is more gen- 
eral than the Lagrangian, for purposes of the calculation 
we use Eq.fJlJ explicitly. Below we summarize the details 
of the calculation. 



A. Calculation 

To compute the process ip —> ZZ + Zj + 77 — > Al we 
include the diagrams shown in Fig. [2] and parametrize the 
scalar coupling to gauge bosons as in Eq. [TJ The total 
amplitude can be written as, 

M = M zz + Mzy + M-fZ + A4 77 (3) 
which upon squaring gives, 

\M\ 2 = \Mzz? + \M Zl \ 2 + \M. fZ \ 2 + \M 7 y\ 2 
+2Re[MzzM*z 1 +MzzM* lZ +MzzM* 11 (4) 

M 11 M*z 1 + M 11 M* 1 z+Mz~ 1 M* 1 z) ■ 
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all couplings at every vertex encountered when tracing over 
the Dirac strings as distinct as well as all Breit-Wigner fac- 
tors so for any amplitude squared term there can in principal 
be four different vector bosons as intermediate states. In the 
case of the photon we have of course g R = g\ = — e em and 
m 7 = T 7 = 0. Since at this stage the various couplings and 
masses are completely general, Eq. ([6jl applies to any process 
where a scalar decays to two spin-1 vector bosons which then 
decay to massless fermions through 'Z-like' couplings. 

Expanding out the terms in Eq. (J6j) we can write the ampli- 
tude squared as, 



.vi,.vr = C + ±L+± + C+~L + ~ + 



ijij ijij ijij ijij ~ ^-f »JV »j» 



(8) 



FIG. 2. Feynman diagrams contributing to tp — > 2i\2(.2- The 
arrows are to indicate the direction of momentum flow. 



An explicit calculation of all of these terms is overwhelm- 
ing, but things can be simplified greatly by taking the fi- 
nal state leptons to be massless. In this case, the momen- 
tum dependent terms in the Z boson propagator numer- 
ators do not contribute. This leads to the propagators of 
both Z and 7 to have the same Lorentz structure, namely 
the Minkowski metric g^ v . This implies that all of these 
terms have the same general Lorentz structure. The only 
difference from these terms comes from Breit-wigner fac- 
tors in the propagators as well as in the couplings of the 
vector bosons, some of which are zero thus 'turning off' 
the contributions from their corresponding Lorentz struc- 
ture. To see this, let us consider the amplitude for any 
combination of intermediate Z and 7 shown in Fig[2j 



Ma = uj (i-y J (g 3 2R P R + gl L Ph) ) v 3 



-ig v 



M? - m? + imjT, 



lJ \M1 -m? + irrnTi 



u i (ij^igiRPn + 9ilPl)) Vi (5) 



where i,j label Z or 7 while 1 and 2 label the final state lep- 
tons and can in principal be e or /x. In the 4e and 4/i case one 
must also include the interference between identical particles, 
but we do not address that issue here. Upon squaring the am- 
plitude and summing over final state lepton polarizations we 
can obtain a general amplitude squared which encompasses 
any of terms in Eq.(|4| and is given by, 



where 



MuMfj = (DuDtjDliDZ-y 1 
9^Vd ' 9i^){g^,%'jj9^) T tj V jj' y 



Vi = (gl R gl R + gl L glL)Tr(Ar far 

+{9\r9\r- fliiflk)Tr(^i 7 * AtV) 
D u = M\ -ml + iTirm 



(6) 



(7) 



The g R ^ are at this point gen- 



and T^ are given in Eq.jfJ). 
eral left and right handed couplings of a 'Z-like' spin-1 vector 
boson to a pair of fermions. The bars are to indicate that 
the corresponding index belongs to the conjugated amplitude 
and are distinct indices from the un-bared ones. We treat 



where a,b — (+, — ) with a and b corresponding to the fermion 
pairs labeled 1 and 2 respectively and 



(D li D 2j Dl-,D; ] ) 



(9) 



Lf^ = {9^T?ig M ){g VJ T2Zg^)T%T^ 



The Ti± are the Dirac traces found in Eq. Q and ± indicates 
whether the trace ends with a 7 s (— ) or not (+). The full 
amplitude squared can then be built out of the object^] in 
Eq.M, 

{M ii Mf 3 ) ab =C$ 3 L$ s . (10) 

Since all of the angular information is contained in the L°^j 
we can take advantage of the simple nature of these terms to 
perform the desired integration before summing over ± and 
the various vector boson intermediate states, after which ana- 
lytic integration becomes unmanageable. Expressions for the 
L?%? are obtained in terms of invariant dot products and CM 
variables. These objects can be used to build the differential 
cross section of any scalar decay to four massless fermions via 
two spin-1 vector bosons. From these one can also reproduce 
analytic results for other processes such as the semi-leptonic 
decay of the Higgs to lujj [26] . 

The final fully differential decay width can now be written 

as, 



dM'fdMXdil 



= n 



\ ijij ijij 



(11) 



where dSl = dc<s>dce 2 dco 1 dQd<&\ (eg = cos 8) and fl^ is the 
final state lepton four body phase space derived following |27| 
and given by, 



n« = (^) 2 ( 1 - 



2ty' k 32tt 2 ' v 32tiV 



1 + 



(M? - Mi) 2 



2(M? +M|)\i/2 



We can obtain the differential mass spectrurrj^] via 



ab__' 
'ijij- 



(12) 



(13) 



1 Expressions for the various coefficients and Lorentz structure can 
be obtained by emailing the corresponding author. 

2 We give an analytic expression for a particular hypothesis in the 
Appendix. 
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We note that we perform the sum over vector bosons before 
the sum over ± which allows for greater simplification of the 
expressions. We can now go on to examine the differential 
mass spectrum for different signal hypothesis. In the Ap- 
pendix we show various singly and doubly differential spectra 
for a number of signal hypotheses. We also give in Eq. (|31[) 
of the Appendix, an explicit expression for the doubly differ- 
ential mass spectrum of a scalar with SM-like ZZ couplings 
and both CP even and CP odd Z'y couplings including all 
interference effects. 



B. The Differential Mass Spectra 

In this section we examine the singly differential mass spec- 
tra for various signal hypotheses and give a feel for how Mi 
and M 2 might be able to distinguish between the different o 
erators in Eq.|2j|. Explicitly we consider the following caser* 

• 1: SM including Z r y and 77 (Aizz = 2,A 2Z ry 



0.014, A 2l 



-0.008 



2: SM coupling to ZZ plus enhanced Z-y and 77 
(Axzz = 2, A 2Zl = 3 * 0.014, A 2ll = -1.3 * 0.008) 

3: SM coupling to ZZ plus CP odd couplings to 77 
and Z7 (A 1ZZ = 2,A 3Z -y = 0.01, As-y-y = 0.01) 

4: CP odd/even mixed coupling to ZZ only (A\ ZZ = 
2, Az ZZ 



0.1) 



0.1,A 2ZZ — l,A 2Zl 
l,A 3Z ~, = 0.01, A 3ll 



• 5: General Scalar (A\ ZZ 
0.01, A 2rr = 0.01, A 3ZZ -- 
0.01) 

where we also show the values for the couplings chosen in 
Eq.|l]). Couplings whose values are not shown in a given 
hypotheses are taken to be zero and we take all values at 
-y/s = mh — 125 GeV. Note that all of these couplings can be 
interpreted in terms of the couplings in Eq. |2j| if we assume 
only up to dimension 5 operators contribute. 

We obtain the differential mass spectra via Eq.(13l, fol- 



lowed by integration over Mi or M 2 , and compare them for 
different hypotheses. These are shown in Fig.[3]for two ranges. 
The first range we take 4 GeV< Mi, a < 120 GeV treating M x 
and M 2 symmetrically shown in the top plot. In this case we 
only show the Mi distribution since it is identical to the M 2 
distribution and only show the lower mass region above which 
the different cases are very similar. 

We also consider the more 'experimental' cut requiring a 
wide window around the Z boson mass 40 GeV < Mi < 
120 GeV and 4 GeV < M 2 < 120 GeV for the 'off-shell' vector 
boson. In this case the Mi distribution is indistinguishable 
for the separate cases so we only show the M 2 distribution. 
One can see, that in particular in the low mass region, these 
variables can be highly discriminating between the different 
cases. We point out also that, if values of M 2 < 10 GeV can 
be probed, the requirement of a window around a Z boson 
may decrease sensitivity to certain hypotheses which have a 
sizable 77 or Z7 component such as hypothesis 5. 



3 We have validated these cases with FeynRules/CalcHEP [161128] 
and the Monte Carlo generator introduced in [3j. 

4 Values obtained from 25 after translating to our parametrazia- 
tion. 
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FIG. 3. In the top plot we take 4 GeV < Mi, 2 < 120 GeV 
while in the bottom plot we take the range 40 GeV < Mi < 
120 GeV, 4 GeV < M 2 < 120 GeV at </s = m h = 125 
GeV when integrating over phase space. The SM is shown 
in blue, but is essentially indistinguishable from hypothesis 3 
(see text). 



Our lower bound on M 2 is chosen to be 4 GeV since lower 
values of M 2 runs the risk of contamination from J/ip states 
whose mass is ~ 3 GeV. We emphasize that experimental 
analyses should be made to push down as far possible since 
as can be seen in Fig. [3] one needs to be able to probe M 2 
below ~ 10 GeV in order to discriminate between a SM scalar 
(hypothesis 1) and one with enhanced Z7 and 77 couplings 
(hypothesis 2) for example. Though current experimental sig- 
nal searches in the golden do not yet consider such low values 
for M 2 , it seems feasible to push the M 2 cut down further as 
was done in the CMS observation of the Z — > £ + £~£ + £~ pro- 
cess [25]. We therefore include this highly interesting region 
here and hope that it may motivate efforts to push the M 2 
reach lower. We leave a complete analysis including detector 
effects to an ongoing study. 



IV. BACKGROUND 

The dominant irreducible background to the golden chan- 
nel comes from qq annihilation into gauge bosons. At energies 
~ 125 GeV the dominant contribution comes from t-channel 
Z7 production. However, as we will see contributions from s- 
channel Z — > 4£ diagrams can effect the angular distributions 
such as the distribution of the angle between the decay planes 
$ defined in SeclIT] Furthermore, we include the ZZ and 77 
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contributions since in principal these are always present and 
may have observable interference effects due to the fact that 
they add at the amplitude level when decaying to charged lep- 
tons. In addition, the inclusion of these contributions allows 
for considering a much larger energy range in one fully dif- 
ferential cross section than can be considered when including 
only the t and u channel contributions. Of course NLO ef- 
fects, including the gg initiated process [30H32J will contribute 
as well, but these are expected to be small and mainly only 
effect the 'input' invariant mass (and overall normalization) 
for the fully differential cross sections. We will examine this 
point below. 

It should also be noted that ideally one would like to include 
the 4e and 4fi final states which in some kinematic regimes can 
have non-negligible contributions from interference between 
final state particles |18| . The inclusion of this channel would 
allow for greater sensitivity for the same amount of luminosity. 
However, because of the interference between identical final 
states in this case, the Lorentz structure becomes severely 
more complicated and we thus leave this calculation for future 
work. 



q fi q 




FIG. 4. Feynman diagrams contributing to qq — > ZZ — > 2e2\i 
and qq — > Z — > 2e2/i. The arrows are to indicate the direction 
of momentum flow. 



A. Calculation 

The background calculation is much more involved than 
the signal calculation due to a higher number of Feynman 
diagrams in addition to a more complicated Lorentz structure. 
As in the signal case the amplitude can be written as, 

M = Mzz + Mz-y + M~,z +X 77 . (14) 

Now however, each of these amplitudes breaks down into six 
'sub-amplitudes'. To see this, let us first consider the ZZ me- 
diated decays. There are three diagrams which contribute to 
the 2e2/i process shown in Fig. [4] First there is the t-channcl 
contribution shown in the bottom diagram. This contribution 
(and its u-channel counterpart) has been computed previously 
for both on-shell |21j and off-shell [5] Z bosons. The second 
contribution comes from resonant 2e2p production proceed- 
ing through the top two diagrams. Each of these diagrams 
also has a corresponding 'crossed' diagram taking into account 
the other possibility for attaching the vector boson lines. This 
gives six diagrams for the ZZ contribution to the golden chan- 
nel. Similarly, there are six more for the 77 contribution plus 
six for Z~y and six for 7Z giving a total of twenty four dia- 
grams. At first this many diagrams can seem intractable, but 
as we will see, when organized in a proper manner the calcu- 
lation is relatively straightforward with the help of Tracer |33| 
and Mathematica [34] to perform the Lorentz contraction and 
symbolic manipulation. 

To begin we first note that the six diagrams can be 'twisted' 
and arranged into the form found in Fig. [5] where we now al- 
low the vector bosons to take on any Z or 7, but once chosen 
are treated as fixed. We use the conventions indicated in 
the diagrams and in particular refer to the diagrams (a), (c), 
and (e) as 't-channel' type diagrams and (b), (d), and (/) as 
'u-channel'. This is not to be confused with the typical vo- 
cabulary for this process which refers to diagrams (a) and (b) 
as t and u channel and diagrams (c) — (/) as s-channel. We 
find this re-naming convenient for organizing and reducing 
the many terms which need to be computed for the differen- 
tial cross section. The Lorentz structure for all of these am- 
plitudes is clearly the same. One needs only to keep proper 




(c) 3 1 (d) 




(e) ' 2 (f) 

FIG. 5. Feynman diagrams contributing to qq ViVj — > 
2£i 2£ 2 and qq ->• Vi -> 2£i2£ 2 . Note that diagrams (c) - (/) 
are in fact s-channel diagrams so the fermions labeled by 1 
and 2 are not to be confused as being in the initial state. 
This is taken into account in how the various momenta are 
assigned as indicated by the arrows. 



track of how the various momentum are routed through each 
diagram. We can see this by considering the amplitude ex- 
plicitly. Using the massless initial quark and final state lepton 
approximation we can write any of the amplitudes in Fig. [5] 
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Using the conventions just described, we can write the 
product of any two amplitudes within a row as 



9 3 zr 



M Xij = (n 

«x (i^(9xRPR + 9x L PL)) ( 



-m 2 +im,T. 



^ (Wx^ + <?x^))«x 



(15) 



where we label the amplitude by the 'long' dirac string, in this 
case X. The labels X/Y/Z — 1,2, q where 1,2 are for final 
state lepton pairs while q is for the initial state quarks. The 
i,j — Z,"y label the vector bosons and n — t,u labels the t 
and u-channel diagrams in our new vocabulary. The invariant 
masses are defined as My = (py +p Y ) 2 - The internal fermion 
momentum are given in terms of external momentum by, 

kqt =Pq~ (Pi + Pl), k qu =Pq- (P2 +P2) 

kit = -Pi - (j)2 +P2), klu = (p q + Pq) - pi 

k 2t = (Pq + Pq) - P2, k 2u = -P2 - (Pi + Pi) ■ (16) 

Note that the invariant masses My and Mz do not necessar- 
ily correspond to the invariant mass formed by the final state 
lepton pairs, as they do in the signal case and in previous 
analytic calculations of the golden channel background which 
neglect the s-channel diagrams. Now with the inclusion of the 
resonant four lepton processes in (c) — (/) we have for these 
diagrams M Y / Z ~ Mq = s which is equal to the invariant 
mass of the four lepton system. To obtain any of the phys- 
ical am plitudes one simply assigns the appropriate labels to 
Eq.(15l as well as the appropriate momenta. Thus for exam- 
ple for diagram (c) we have X — > 1, Y — ¥ q, Z — ¥ 2, and 
n — > t. To switch from t-channel type to u-channel diagrams 
(staying in the same row in Fig[5| one simply takes t — > u 
and 7 CT o 7 7 . Of course at this stage all these labels are 



arbitrary meaning that the amplitude in Eq.(15l applies to 
any process with this topology and Lorentz structure. Note 
that for the Z propagators we drop the momentum dependent 
terms since they do not contribute in the massless lepton ap- 
proximation. As mentioned in the signal case, for the photon 
9r = 9l an d m i = T 7 = 0, but for now we take the couplings 
and propagators as general. 

As in the case of the signal, the next step is to find a gen- 
eralized amplitude squared for any two of the six diagrams. 
Although there are in principal thirty six terms when squar- 
ing the amplitudes, these organize themselves into only two 
distinct types of Lorentz structure. The first type is found 
when multiplying any two diagrams in the same row of Fig. [5] 
This is the Lorentz structure found in our previous calcula- 
tions of the ZZ contribution in which only diagrams M q t and 
A4 qu are included (the top row). The square of Mu + A4i u 
and Mit + M.2u (second and third rows) will also exhibit 
this Lorentz structure. The second type of Lorentz structure 
is obtained when taking the product of any two diagrams 
in different rows. In the conventional language, interference 
between the first and second row or first and third row cor- 
responds to interference between the t-channel di-boson pro- 
duction amplitudes and the s-channel diagrams. Interference 
between the second and third row corresponds to interference 
between the two types of s-channel diagrams. We first dis- 
cuss the 'squared' terms where the amplitudes are contained 
within the same row before examining the interference terms 
between rows. 



M xij M% = (Dyi Dzj D Y iD* Z j ) ~ 1 
(a T' T ' T -a--)(n T 77 -o-- lT" M '" y 



(17) 



where the Tyii an< A Dyi are defined similarly to those in 
Eq.Q and the long Dirac string is given by, 



^xijJJnm ~ (9xr9xr9xr9xr + 9xl9xl9xl9xl) 
(dx r9x r9x r9x r ~ 9xl9xl9xl9xl) 



(18) 



Expanding out the terms in Eq.( |17[ l we can organize in a 
manner similar to Eq.|8| and write the amplitude squared as, 



\An \ a Tn * v~~* nabc j abc 

MxijMxij — l^^XXijij^XXnn 



(19) 



where again a,b,c — (+, — ) in the order X, Y, Z and 

(9xr9xr9xr9xr + 9xl9xl9xl9xl) 
\9yr3yr ± 9yl9yl)(9zr9zr ± 9zl9zl) (20) 

-^xxnm = (9^T§±gv<T){g^Tz\g^)T x p f m± . 

The T Y a z ± are the Dirac traces found in Eq. uTp while the 



-Txrami are those found in (18 1. Again ± indicates whether 



the trace ends with a 7 5 (— ) or not (+). We note that unlike 
in the signal case, when organized in this way (essentially by 
powers of 7 5 ) the gauge structure completely factors from the 
Lorentz structure. This allows us to sum over all possible 
intermediate vector bosons at this stage to write, 



M X MT = T,_M xij M^ 3 



= E 



Epabc r abc 

^XXijij^XX; 



(21) 



E(Ec ai " 



XXijijJ^XXnm 



EC 

abc 



abc T abc 
XX^XXnm 



This simplifies things greatly and in particular the objects 
Cxx now contain all of the information regarding the interme- 
diate vector bosons including the interference effects between 
the different processes. These will serve as overall coefficients 
for the various Lorentz structure pieces. 

We are now in a position to examine the interference terms. 
Let us take the product of any two diagrams not in the same 
row. One can show explicitly, 



M\,\r;:, = {D Yi D Zj _D ZJ D XJ r^ 

(9»<TT x ^- jn g^){g vl T Y l il - jm gp,v)T2 :j - i 



(22) 



where the T2° are as before and the new Dirac strings are 
given by, 



'Xijj 



n — (9xr9xr9xr + 9xl9xl9xl) 
Tr(i, xl v Jtx n ^ i>x^) + 

(9xr9xr9xr ~~ 9xl9xl9xl) 
•2>(i>x 7 " ^„7" t>xl'f) • 



(23) 



The distinct Lorentz structure found here as compared to 



that found in Eq.( 171 is due to the different path taken when 



tracing over the fermonic strings. 



Again we expand out the terms in Eq.(22l to obtain 



A A n A/I m * nabc r abc 

JVl Xi] JVl YTj — 2^ L XYiji] Ij XYr, 
abc 



(24) 



where, 



C 



±±± 



V i - - (D Yi D zj Dl i D* x j)-\gi R g^ R ±g^ L ) 
xr9xr9xr i 9xl9xl9xl)(9yr9yr9yr ^ 9yl9yl9yl) 



r±±± —In T"^ n—\( n f^ 1 " 7 n- - ^T 1 " 
^XYnm ~ \9fJ-aJ-Xn 9^1 AH"1 A Y m Sc" I A Z ■ 



(25) 



and the T x ^ are the traces found in Eq.(23l. As mentioned 



above, since the gauge and Lorentz structure factor com- 
pletely we are free to perform the sum over the intermediate 
vector bosons at this stage once again to obtain the various 
Lorentz structure coefficients, 



M\MT = V An,;.Vi;< : ; 



E, r>abc \ r abc 
\ Z_, L-XYifij J^XYnm 



a be 



abc t abc 
XY^XYnm 



(26) 



Thus again all of the information concerning the intermedi- 
ate vector bosons is contained in C XY - We now have all of 
the piece^] necessary to build the total amplitude squared of 
the diagrams in Fig. [5] including all contributions from the 
intermediate vector bosons. Explicitly we have, 



E ((c$ 



be 
11 



\M q - 

Lobe 
aant 



Mi +M 2 \ 
+ C11 



abc t abc . pabc r abc 
-^llnm~T^22 lj 22nm 



= E 

abc nm 

i r) Ty ( pabc r abc , pabc r abc ■ pabc r abc 
-\-Zrte\Uqi iiqljim t tl2 -^12nm ~r ^2q ^2qnm 



)) 



(27) 



where the sum over intermediate vector bosons has been al- 
ready implicitly performed and the sum over n, m which in- 
cludes the t and u channel contributions is shown explicitly 
(note that this also factors from the vector boson sum). The 
Cxy coefficients are in general complex due to the factor of 
i multiplying the decay width in the massive vector boson 
propogators. The Lorentz structure is either purely real or 
purely imaginary depending on whether the term contains an 
even or odd number of traces ending in 7 s . These traces give 
an overall factor of i (and an epsilon tensor). Thus if L^ Ynm 
contains an even number of these traces, then it is purely real 
and if it contains an odd number it is purely imaginary. The 
squared Lorentz structure Z/xxnm however is strictly real as 
are the squared coefficients Cxx- Taking this into account, 
we can write for Eq.(27l the final amplitude squared as, 



\Mu\ = \M q + Mi +M2 

even 

■ J qqR^- J qqR 



\ ^ / pabc j abc . pabc j abc 
/ ^ \ ^qqR^qqR ~T ^WRL'WR 



abc 
even 



pabc j abc 
^22Rp22R 



pabc j abc \ 
^2qR Li 2qR J 



. X ^ / pabc t- abc . pabc j abc 
T / J \ ^qlR^qlR + ^12R^12E 

abc 

odd 

Ef pabc j abc . pabc j abc . pabc j abc \ 
{^qll^qll +^127-^12/ ~T ^2ql ^2ql j 



(28) 



where we have now performed the sum over t and u channel 
diagrams and Cf YRJ = Cf Y ± C"y respectively. We have 
also implicitly included a factor of 1/4 from averaging over 
initial state quark spins and a color factor of 1/3. The sums 

labeled even = (+ + +,H , — I — , h) indicate terms 

with even powers of 7 s and those with odd = ( h +, H 

+ , + H — , ) indicate terms with odd powers of 7 . Note 

that since the photon has vector like couplings where gL = 9r 



are zero for the 77 

' + 



all coefficients Cxy with a, b, or c 

intermediate state. Thus 77 only contributes to the C XY 
coefficients (including of course when X = Y). 

Previous calculations of the golden channel background, 
which include only the di-boson productio n pr ocess, are con- 



tained within the first term C qq c Lf q c of Eq.(|28|). All the other 



terms arise from the resonant four lepton production process 
and the interference between it and the di-boson production 
process. Note that Eq.(|28[) is also more general than for just 
the golden channel. In principal this expression holds for 
any process with the same topology and 'Z-like' couplings 
to fermions. Since we have built the expression out of a gen- 
eralized Lorentz structure with coefficients, it can easily be 
adapted to consider new physics contributions which may en- 
ter with the same topology and alter some of the coefficients 
by an observable amount. Thus one can imagine performing 
stringent tests of the SM using this parametrization to ex- 
tract the various Lorentz structure coefficients. We leave an 
investigation of this to future work. 

The final fully differential cross sectioij^jis again obtained 
by combining the amplitude squared with the invariant four 
body phase space (see Eq.|12|), 



= U4AM4 



dM^dM%d(l 

The differential mass spectrum^] is obtained again via, 



(29) 



(30) 



We now examine how the various components of the back- 
ground contribute to the differential mass spectrum. 



B. The Differential Spectra 

In this section we examine how the individual components 
of the background contribute to the invariant mass spectrum 
of the four lepton system. In addition we also study how in- 
cluding parton distribution functions (pdfs) and NLO correc- 
tions change the differential spectra by comparing normalized 
projections obtained from our analytic expression to Monte 
Carlo generated by POWHEG [35H38] and Madgraph [17]. 

We first separate the background into its various compo- 
nents which we define as the following, 

• A: s-channel 2e2u process 

• B: t+ u- channel 77 

• C: t + u-channel ZZ 



Expressions for the various coefficients and Lorentz structure can 
be obtained by emailing the corresponding author. 



6 This expression has been validated with the Madgraph matrix 
element squared. 

7 An analytic expression for the dominant component to the back- 
ground is given in Eq.l33|l of the Appendix. 
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• D: t + -u-channel Z*y 

• E: t + it-channel interference only 

• F: ZZ + Z*y + 77 s/t-channel interference only 

where now s, t, and u are used in the usual sense and the res- 
onant s-channel 2e2fi process can proceed through any com- 
bination of Z and 7. We first consider the relative fractions 
of these components as a function of the invariant mass of the 
four lepton system for the range 100 — 600 GeV in Fig. [6] The 
dotted lines indicate when a contribution is negative, which 
of course only occurs for interference terms in certain energy 
ranges when the interference is destructive. The solid black 
line at constant value of 1 is the total partonic level qq — > 2e2/x 
(q = u, d) background including all interference and all inter- 
mediate vector bosons. From Fig. [6] one can see how the 
relative contributions coming from the different components 
change as a function of energy. 

Component C (the ZZ t + u channel) is the only piece of 
the background to have been previously calculated analyti- 
cally [5T] . This makes up the dominant contribution above 
the ZZ threshold, but is negligible from 110 GeV < y/s < 140 
GeV and in fact is even smaller than the interference terms. 
We also plot the spectrum if one requires a window around 
the Z boson mass in the bottom plot of Fig. [6] The dominant 
component near the resonance mass of 125 GeV is D regard- 
less of the window on the Z mass. Except for component F, 
one can see that the relative fractions are fairly insensitive to 
the Z window requirement except in the range ~ 100 — 110 
GeV. 

The flexibility of the analytic expressions also allow us to 
easily isolate the contribution coming from interference terms. 
Component E for example is due to the interference between 
the intermediate gauge bosons in the t + u channel and is 
destructive over the entire range regardless of the Z window. 
The interference between the resonant s-channel and the t- 
channel pair production processes is shown in F and switches 
between constructive and destructive if one requires a window 
around the Z, but otherwise is constructive. Though these 
components are small it is possible for them to have subtle 
effects on the angular distributions such as in the modulation 
of the azimuthal angle <f> (See Fig. 13 1 and may be partic- 
ularly interesting to study in the range 100 GeV < y/s < 
110 GeV. The expressions for most of the components them- 
selves are too cumbersome to write here, but in the Appendix 
we give the expression for the doubly differential (Mi,M2) 
mass spectrum of the full t + u (the sum of B-E) component 
which as we can see in Fig. [5] and Fig. [7] provides a very good 
approximation above y/s ~ 110 GeV. 

To examine the effects of NLO contributions and pdfs we 
compare our parton level result for qq — > 2e2/i (g = u, d) to 
Monte Carlo data generated by the NLO POWEG and LO 
Madgraph codes which include pdfs [39]. For this we define 
our phase space as 40 GeV < Mi < 120 GeV and 10 GeV 
< M 2 < 120 GeV for the energy range 110 GeV< y/s < 
140 GeV. We also plot the t + u component only (defined 
as the sum of B-E) to examine what affects neglecting the 
resonant 2e2/i process has. 

In Fig. |7|9| we show the kinematic distributions where it 
can be seen that NLO and pdf contributions affect the nor- 
malized spectra negligibly. In addition we can see that ne- 
glecting the resonant process also has little effect on all the 
kinematic variables except $, where it affects the modula- 
tion and in the forward regions of cos Si. As we will see in 
the Appendix, the modulation is due almost entirely to the 




200 300 
is (GeV) 



600 




) 200 300 400 500 600 

is (GeV) 

FIG. 6. The four lepton system invariant mass spectrum 
(without pdfs) for the various components defined in the text. 
The dotted lines indicate when the interference between com- 
ponents is destructive, thus giving a negative contribution. 
In the top plot we take the ranges 4 GeV< Mi )2 < 120 GeV 
while in the bottom plot we take the range 40 GeV< Mi < 
120 GeV and 10 GeV< M 2 < 120 GeV while taking Mi > 
Mi. 



resonant process. These distributions simply reflect the fact 
that the various kinematic distributions are not highly corre- 
lated with y/s allowing us to take y/S essentially as an input 
from the pdfs. To build a complete hadronic differential cross 
section one could convolve the y/s spectrum obtained from 
Madgraph or POWHEG with the partonic differential cross 
section obtained analytically. This of course is what would 
be done for an LHC analysis, but we do not do that here and 
instead simply integrate our partonic differential cross section 
over y/I. 

From Figs. [6] and [7] we expect the doubly differential spec- 
trum obtained from the t + u component only to be a good 
approximation which could be useful for a simplified analysis. 
We give an explicit expression for this component in Eq. ( 33 1 
of the Appendix. Though it does not use all of the kinematic 
variables, it should still have strong discriminating power and 
can be used with the methods proposed in [7j to form a pow- 
erful simplified study. 
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FIG. 7. Comparing the LO and NLO results for the Mi and 
M2 invariant mass spectra for the ranges 40 GeV < Mi < 
120 GeV and 10 GeV < M 2 < 120 GeV. We take the range of 
the four lepton system invariant mass to be 110 < ^/s < 140 
GeV. 



V. CONCLUSIONS AND OUTLOOK 

We have calculated and presented analytic fully general 
differential cross sections for the golden channel signal and 
background in the 2e2/i final state including all intermediate 
vector bosons and interference effects. We have presented var- 
ious singly and doubly differential spectra and examined how 
the different interference effects manifest themselves in these 
distributions and in correlations between the different kine- 
matic variables. We have also emphasized the need to push 
the 'off-shell' invariant mass (M2) reach as low as possible as 
well as relaxing the 'Z- window' to maximize the discriminat- 
ing power when testing different signal hypothesis. We have 
shown that the expressions can aid in distinguishing between 
different signal hypotheses and because signal and background 
are provided, both can be included into one likelihood, as 
should be done when performing simple hypothesis tests of 
different signals. These expressions can be implemented into 
an MEM analysis to perform detailed studies of the spin and 
CP properties of any scalar resonance which has been or may 
be discovered at the LHC. 
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FIG. 8. Comparing the LO and NLO results for the polar 
angles cosO, cos#i, cos #2 for the ranges 40 GeV < Mi < 
120 GeV and 10 GeV < M 2 < 120 GeV. We take the range of 
the four lepton system invariant mass to be 110 < -y/s < 140 
GeV. 
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FIG. 9. Comparing the LO and NLO results for the azimuthal 
angles $ and 3>i for the ranges 40 GeV < Ah < 120 GeV and 
10 GeV < M 2 < 120 GeV. We take the range of the four 
lepton system invariant mass to be 110 < ^/s < 140 GeV. 



VI. APPENDIX 

The general scalar and background differential spectra are 
too cumbersome to write in one pag^] for most of the differ- 
ent components. We give a couple of the simplest ones here, 
but that are not found in literature. We also examine how 
the different signal hypotheses and background components 
contribute to the various kinematic distributions. We show 
a multitude of singly and doubly differential distributions for 
both signal and background. Of course none of these plots can 
show the discriminating power of the fully differential cross 
section, but one can visually get a sense for the discriminat- 
ing power of these kinematic variables. Detector effects will 
also shape these distributions and deserve careful study, but 
it is clear that the golden channel is a powerful probe of the 
underlying physics. 



A. Analytic Expressions 

We give here a pair of analytic expressions for the differ- 
ential mass spectra for one of the signal and one of the back- 
ground components which are simple enough to fit on one 
page. Although not as powerful as using the fully differential 
cross section, with just these two relatively simple expressions 
one can perform robust analyses of the newly discovered scalar 
and its coupling to neutral gauge bosons as suggested in [7]. 
For the signal we give the ip — > ZZ + Z*y — > 2e2/i differential 
mass spectrum including interference. For the ZZ coupling 
we take only the 'SM-like' coupling A\zz to be non-zero. For 
the Z7 coupling we allow for both A 2 z-y and Azz-i to be non- 
zero, thus allowing for CP violation. Using Eq.( 13 1 we obtain, 



dTsM+z-j 
dMldMl 



Mf + (Mf - s) 2 - 2M 1 2 (Mf + s)(6A 1Z zA 2 zAe l (g L + g R )(gl + g 2 R )M 1 M 2 m 2 z (2M 2 Ml - (M 2 + M 2 2 )m 2 )(-1 + /3 2 ) 
(1 + + Pi) + A 2 lzz (g 2 L + g\fM\Mlm\^\ - J\^[\ - /3 2 2 (3 + Ifafo - 2ft 2 + /3 2 (-2 + 3/3 2 2 )) + 2e 2 ^/l - ftyfl - /3 2 

(2A 2 ZA (/?i + /3 2 ) 2 + A 2 2ZA (3 - ft 2 + 4Pifc + (-1 + 3/3 1 2 )/3 2 2 ))(2 ffiffii M 2 M|((Af 1 - m z )(M 2 - m,)(Afi + m z )(M 2 + m z ) + m 2 z Y 2 z ) 
+gl(M^m 2 z (m 2 z + T 2 ) + M 2 M 2 m 2 ( — 3Mf + m 2 + T 2 ) + M^Mo 1 - 3M 2 2 m 2 + m\ + m 2 F 2 )) + <?fl(M 2 4 m 2 (m 2 + T 2 ) 
+M 1 2 M 2 2 m 2 (-3M 2 2 + m 2 + T 2 ) + JW?(3Af£ - 3M 2 2 m 2 + m\ + m 2 z T 2 z ))] 



'(46087r 4 s 2 ^(l - /3 2 ) 3/2 (l - /3 2 ) 3/2 ((M 2 



m 2 z ) 2 + m 2 F 2 )((Mf - m 2 ) 2 + m 2 r 2 )j 
I 



(31) 



where we define, 
£1,2 = 



4M 2 



^ (l±(M?-Mf)/ s y 



(32) 



This expression is frame invariant and can accommodate a 
Higgs-like particle with SM couplings to ZZ, but with perhaps 



new physics contributions through its couplings to Z"/. The 
ei are the photon couplings to charged leptons while gL,R 
are the leptonic Z couplings. Mi and M 2 are the final state 
lepton pair invariant masses while mz is the mass of the Z 
boson and -^/s is the four lepton system invariant mass. The 
doubly differential mass spectrum for the full t + u component 
of the background (sum of components B-E) can be obtained 
analytically via Eq. |30| to give, 



The distributions will be made public in the near future, but can 
be obtained from the corresponding author in the meantime. 
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da 



BG 
't + u 

dMfdMi 



foU + 9ta)(d + 9H?M£M* + 2e,e 9 (<4 + g 3 qR )(g L + g R )(gl + g R )MlMl(2MlMl - (M( + Ml)m 2 z ) 
+ 8efe q ((M 2 - mlf + m 2 F 2 )((M 2 2 - m 2 z ) 2 + m 2 r 2 ) + Aefe q (g qL + g qR )(g L + g R )((M 2 - m z )(M 2 + m,)(-Af? + m 2 z ) 
{-2M?Mi + (M 2 + M|)ro 2 ) + m 2 z {Mt + M% - (M? + M 2 2 )m 2 )r 2 ) + 2e 2 e 2 q {g 2 qL + g 2 qR )(4g L g R M?M!(M 1 - m z )(M 2 - m,) 
(Mi + m z )(M 2 + m.) + ,g£((-2M 1 2 M 2 2 + (A/ 2 + M|)m 2 ) 2 + (M 4 + M 2 4 )m 2 r 2 ) + ff |((-2M 1 2 M 2 2 + (M 2 + M 2 2 ) 

m 2 ) 2 + (M 4 + M 2 4 )m 2 r 2 )))(4(M 1 2 + M 2 2 - s^A/ 4 + (Mf - s) 2 - 2A/ 2 (Mf + s) - ((M? + M 2 2 ) 2 + s 2 ) 

(log [(A/? + M| - s + sJm* + (Ml - s) 2 - 2M 2 (M| + s)) 2 ] - 21og [- Mj 2 - M| + s + sj + (Mf - s) 2 - 2M x 2 (Mf + «)]))) 
/ (27648 M 2 M|7r s (Mi + M| - s)s 2 ((A/ 1 2 - m 2 ) 2 + m 2 r 2 )((A/| - m 2 ) 2 + m 2 r 2 )) . (33) 



This expression includes the ZZ, Z^ and 77 contributions 
including all interference and can be combined with pdfs or 
be used for a leptonic initial state. The e q are the photon 
couplings to the initial state fermions while the g qR /L are the 
initial state fermion couplings to Z bosons. Note that these 
expressions have not been normalized and should be thought 
of as at fixed s. 



B. Singly Differential Angular Distributions 



C. Doubly Differential spectra 



In Fig. [I0p4l we show the angular distributions for the 
5 angles (cos O, cos 9\, cos 82, $1, 5?) found in the four lepton 
system and defined in Sec. [Il] We plot the angular distribu- 
tions for signal hypotheses 1-5 defined in Sec. |IIIBl and also 
show the various background components defined in Sec. |IV B] 
For all distributions the phase space is defined as 4 GeV 
< Mi < 120 GeV and 4 GeV < M 2 < 120 GeV with 
Vs = 125 GeV for signal and 110 GeV < %/S < 140 GeV 
for background. 

Since we are considering a spin-0 scalar as our signal, the 
cos and $1 are of course flat, but are still useful for dis- 
criminating between signal and background. A particularly 
interesting variable is the azimuthal angle between the lepton 
decay planes, $. This is especially sensitive to the various 
interference effects as well as the CP properties of the decay- 
ing scalar, as was pointed out in |22) . One can see that the 
different signal hypothesis affect the modulation of $ while 
an extreme case like the CP violating hypothesis 5 can lead 
to a striking signal in the form of an asymmetric modulation 
and phase shift relative to the SM prediction. 

For the background we can see how the various components 
contribute to the different kinematic variables. It is clear that 
the Zj t + u component (D) is the dominant contribution for 
our defined phase space. Note however, that the s-channel 
component (A) also contributes and in particular is the dom- 
inant contribution to the modulation of $. We can also see 
that the resonant process affects cos 9\ and cos 82 , especially 
in the forward regions. It is also interesting to comment that 
the 77 contribution (B) is featureless in all of the distribu- 
tions except for a small upward slope in the extreme forward 
directions of cos©. Note that for the $1 azimuthal angle, the 
modulation is due entirely to the Z"/ t + u component (D). 
Whether these different effects can still be seen once detector 
effects are included requires careful study which we leave for 
future work. 



In Fig. |15|21| we show various combinations of the doubly 
differential spectra for the five signal hypotheses as well as the 
full background. These are primarily for illustration purposes, 
but from these one can get an idea of the correlations between 
the different kinematic variable^] For these plots only the 
five signal hypotheses and the full result for the background 
are shown. For all distributions the phase space is defined as 
4 GeV < Mi < 120 GeV and 4 GeV < Af 2 < 120 GeV with 
VS = 125 GeV for signal and 110 GeV < V§ < 140 GeV for 
background. 



9 We do not show all possible combinations, but any not shown 
here can be obtained by emailing the corresponding author 
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FIG. 10. On the left hand side we have plotted the cos O angular distributions for hypotheses 1-5 (hypothesis 1 = SM) defined 
in Sec |III B] On the right hand side we plot the components A-F of the background defined in Sec |IV B| 
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FIG. 11. On the left hand side we have plotted the cos Si angular distributions for hypotheses 1-5 (hypothesis 1 = SM) defined 
in Sec jHI B] On the right hand side we plot the components A-F of the background defined in Sec jlV B] 
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FIG. 12. On the left hand side we have plotted the cos #2 angular distributions for hypotheses 1-5 (hypothesis 1 = SM) defined 
in Sec jHI B] On the right hand side we plot the components A-F of the background defined in Sec jlV B] 
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FIG. 13. On the left hand side we have plotted the <3? angular distributions for hypotheses 1-5 (hypothesis 1 
Sec jHIBl On the right hand side we plot the components A-F of the background defined in Sec |IV B] 
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FIG. 14. On the left hand side we have plotted the $i angular distributions for hypotheses 1-5 (hypothesis 1 
in Sec jHI B] On the right hand side we plot the components A-F of the background defined in Sec jlV B] 
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FIG. 17. The (M2,3>) doubly differential spectrum. The first five distributions are for signal hypotheses 1-5 (hypothesis 1 = 
SM in top left) defined in Sec IIIB while the bottom right plot is for the full background. 
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FIG. 19. The (cos 9i, cos 62) doubly differential spectrum. The first five distributions are for signal hypotheses 1-5 (hypothesis 
1 = SM in top left) defined in Sec IIIB while the bottom right plot is for the full background. 
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FIG. 21. The ($,$1) doubly differential spectrum. The first five distributions are for signal hypotheses 1-5 (hypothesis 1 = 
SM in top left) defined in Sec IIIB while the bottom right plot is for the full background. 
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